La7/8Sri/gMn03: long period orbital order with hole stripes 
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The electronic structure of Lay/gSri/gMnOs was studied by using the rotation-invariant 
LSDA (Local Spin Density Approximation) +U formalism. We found the specific kind of charge and 
orbital order for ferromagnetic insulating solution in 2ax26x4c supercell, where a, b, and c denote 
cubic perovskite axes. The holes are distributed over 3-line stripes in every second ac plane along 
a axis. The stabilization mechanism of the hole stripe is attributed to the strong p-d hybridization. 



75.50.Dd, 71.20.Be, 75.30.Et 

Much attention has been focused on doped perovskite 
manganite systems Lai-^Sr^MnOs, because of its metal- 
insulator transition, colossal magnetoresistance, spin 
structures and properties, and charge and orbital or- 
der. The lightly doped Lay/gSr^gMnOs shows the 
reentrant structural transition of orthorhombic structure 
O* — >0' — >0*, where the transition temperatures are 
160 K and 260 K, respectively pL Further, the sam- 
ple shows the ferromagnetic transition at T c =200 K [Q. 
The ferromagnetic order suppresses drastically the static 
Jahn- Teller distortion. The resistivity in the O' phase 
below T c shows the metallic temperature dependence 
and, therefore, the boundary of the low temperature O* 
and O' phases coincides with the insulator-metal phase 
boundary. The paramagnetic state is insulating. Yamada 
et al. ||] reported the observation of satellite spots in the 
low temperature O* phase, corresponding to 2ax26x4c 
supercell structure jj), and proposed the polaron order 
in this phase. Endoh et al. Q confirmed the super- 
lattice structure in the low temperature O* phase of 
Lao.88Sro.i2MnC>3 but reported no signal for charge order 
in the x-ray study. They also reported the G-type orbital 
order in the low temperature O* phase, deduced from the 
resonant x-ray scattering. Mizokawa et al. 0] found sev- 
eral stable spin, orbital and charge ordered states, using 
the model Hartree-Fock Hamiltonian. They introduced 
small Jahn- Teller distortion for opening the band gap. 
However the most striking feature of La7/ 8 Sr 1 / 8 Mn03 
may be the observation of orbital order and insulating 
property in a ferromagnetic system without any (or very 
small) Jahn- Teller distortion. 

In this Letter, we present the results of ab-initio 
self-consistent calculations of the electronic structure of 
La 7 / 8 Sr 1 /gMn03, obtained by using the LSDA+U for- 
malism B in the framework of linear muffin-tin orbital 
method in the atomic sphere approximation (LMTO- 
ASA) (7), and the experimentally determined crystalline 
data of low temperature O* phase (|]. The Coulomb 
and the exchange parameters are set to be U = 8.0 eV 
and J = 0.88 eV for Mn al- and U = 6.0 eV for O p-shell, 



respectively |J . The inclusion of LSDA+U potential cor- 
rection for O p-shell was found to be crucial in the present 
calculations, since it controls the value of charge transfer 
energy between Mn d- and O p- valence orbitals |h]] . 

In order to understand the nature of the doped ma- 
terial, we have calculated the electronic structure of 
the model ferromagnetic undoped LaMn03 compound 
with the crystal structure parameters corresponding to 
the doped case. It was found that one e g major- 
ity spin electron per Mn site occupies alternatively 
r)5(03z 2 -r 2 + <f>x 2 -y 2 ) and ^(</> 3z 2_ r 2 - (f> x 2 _ y2 ) or- 
bitals 0JT| (Figjl] (i)). The band gap of 0.58 eV is 
opened in the energy spectrum between majority spin 
states (Figj2j), since there is already an extended gap of 
4.42 eV in the minority spin subband in the vicinity of the 
Fermi energy. One can suppose that unusual ferromag- 
netic insulating state of model undoped LaMn03 com- 
pound is stabilized due to fine balance between double- 
and super-exchange interactions. Without orbital order 
the double-exchange mechanism would bring the ferro- 
magnetic metal ground state. The insulating state ap- 
pears due to the orbital polarization in the twofold de- 
generated e g majority spin band occupied by one elec- 
tron. The super-exchange interaction in orbitally or- 
dered state is ferromagnetic due to specific overlapping 

Of ^(</>3z 2 -r 2 + 0x 2 -y 2 ) and -i=(</>3 z 2_ r 2 - 4>x*-y*) OT- 

bitals on neighboring Mn ions on ac plane. That is an 
example of an antiorbital but ferromagnetic order (for re- 
view of mutual correspondence between orbital order and 
preferable magnetic interactions in insulators containing 
orbitally degenerate transition metal ions see, e.g. fbl). 
The alignment of the same e g orbitals along b direction 
(see Fig. [j] (i)) gives rise to the larger band width and 
resultant further stability than the case of the alignment 
of the alternative e g orbitals. 

To investigate the influence of the Sr doping, we have 
calculated the energy spectrum of Lay/gSrjygMnOa for 
the minimal possible unit cell containing eight Mn atoms 
and one hole. It was found to be a ferromagnetic half- 
metal with low density of states at the Fermi level equal 
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FIG. 1. Angular distributions of e 9 electron spin density (ADESD) of Mn atoms and corresponding spatial orbital order in 
the 2a x 26 x 4c unit cell of Lay/gSrx/gMnOa compound for undoped (i), doped metallic (ii), and doped insulating (iii) solutions. 
For case (iii) there are translated atoms from the top to the bottom plane along c axis, and there are no translated atoms along 
a and 6 axes. The unit cells for (i) and (ii) cases are increased up to the same size as for (iii) case to make the comparison 
more convenient. The ADESD is defined as p(9,tf>) = J2 m ,m> Qm,m'Ym(0,<f>)Y m t(0,4>), where Q m , m / = n] n m , - n l m m , is 2x2 
e 9 -orbital occupation matrix obtained in our self-consistent calculation and Y m (6,(f>) are corresponding spherical harmonics. 



to 0.10 states/(eV-atom Mn). One would expect the uni- 
form distribution of hole state over all eight Mn atoms 
for the metallic solution. However it is not the case. The 
Mnl atom (for denotations see Fig. [l] (ii)) shows the max- 
imum loss of its d-electron density (0.06 e) and reduces 
its d-magnetic moment by 0.10 [Lb in comparison with 
undoped state Jl4| . This deficiency of (i-electron density 
looks as the decrease of the e g orbital size in Fig. [l] (ii) . 
Thus all these facts indicate the preferable occupation of 
Mnl site by the hole states. 

The influence of Mnl hole on the nearest sites Mn2 and 
Mn3 (Fig (ii)) is so strong that the changes of Mn2 and 
Mn3 atoms in d-occupancy (0.04 e), d-magnetic moment 
(0.07 hb) and e g orbital density angular redistribution 
allow to consider these sites as those containing the hole 
together with Mnl. At the same time the change of c g 
states of Mnl' atom is mainly in shape. There are no es- 
sential changes of d-occupancy and d-magnetic moment 
in comparison with model undoped case in the electronic 
states of second plane Mn atoms (backplanes in Fig. |l] (i) 
and (ii) look similar). Thus we conclude that one hole 
per eight Mn sites is distributed in such a way that only 
three Mn atoms out of eight change significantly their e g 
electronic shell. 

All these three hole-doped Mn atoms are in the same 
ac plane (frontplane in Fig. m (ii)). In frontplane the 



e g orbitals of the hole-doped Mn atoms overlap strongly 
leading to the formation of a wide band. On this stage - 
for small unit cell restricted by eight Mn atoms - double 
exchange would win super-exchange and the spectrum is 
metallic. It is interesting to note that partial densities 
of states have vanishing values at the Fermi level for Mn 
atoms on the backplane. Thus, the backplane remains 
insulating even for doped metallic solution. 

However, the real La7/gSri/ s Mn03 compound demon- 
strates the insulating properties. We observed the ap- 
pearance of the small (0.03 eV) band gap in the energy 
spectrum of Lay/gSrjygMnOa (Fig. ||) when the unit cell 
was increased up to 2ax26x4c size, where a, 6, and c 
denote cubic perovskite axes. In this ferromagnetic insu- 
lating case there are two holes per sixteen Mn atoms. 

One can see two "impurity" peaks in the density of 
states in the energy region of former band gap for model 
undoped case (Fig. ||) . The low-energy peak has capacity 
precisely equal to 1 and is well separated from the rest 
band. It is mostly originated from Mn3 and Mn4 states 
with the admixture of Mn2 states. The second peak 
is not split enough from the band to be single but its 
capacity reaches 1 in the deep. States of Mnl and Mn2 
contribute mostly to this peak. In these two "impurity" 
peaks, 79 % of hole states are localized on six Mnl-Mn4 
atoms forming a 3-line stripe (Figs. Ill (iii), 0) and their 
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FIG. 2. Total density of states in the vicinity of Fermi en- 
ergy (Ef = 0) for doped insulating La 7 /gSr!/gMn03 (solid 
line; right ordinate) and model undoped LaMnOs (dashed 
line; left ordinate) compounds. 

nearest oxygen neighbors. The distribution of hole states 
between these Mn and O atoms belonging to the stripe 
is 44 % and 35 %, correspondingly. 

The same conclusion - hole states localization on stripe 
- follows from the values of Mn ci-occupancy and d- 
magnetic moment values. The deficiency of d-occupancy 
is 0.07, 0.05, 0.04, 0.04 e and that of d-magnetic moment 
is 0.11, 0.09, 0.07, 0.07 [i B for Mnl, Mn2 Mn3, and Mn4 
atoms, correspondingly; for the rest of Mn atoms these 
decreases are essentially less. 

The comparison of the mutual orientation of e g orbitals 
between undoped and doped solutions also support this 
conclusion. One can see from Table | for frontplane Mn 
atoms that Mn5 and Mn6 atoms contain no holes in their 
majority spin (f ) e g shell. Wave function of e g electron 
of Mn6 atom is the closest to that of undoped case (close 
to equal contributions of 03 Z 2_ r 2 and 4> x 2_ y 2), whereas 
ej electron of Mn5 atom occupies the orbital close to 
(f>3 y 2_ r 2 wave function fll| . Other Mn atoms presented 
in Table | contain the hole because of their reduced ej oc- 
cupancy. The contributions of 4>3z 2 -r 2 an d 4> x 2 - v 2 wave 




FIG. 3. The pattern of hole order for La7/gSri/gMn03 ob- 
tained in our calculations. Gray circles represent Mn sites 
containing the hole states, white circles - Mn 3+ ions. 



functions provide the angular distribution presented in 
Fig. [I] (iii). And as for the model doped metallic case, 
the backplane of Mn atoms is free of holes. The holes are 
localized on Mn atoms in the 3-linc stripes along a axis 
in every second ac plane (Fig. ||) Jl5[ |. 

In LSDA+U formalism, the less occupied is a localized 
orbital the higher is its energy position Q. The differ- 
ence in e g orbital occupancy of the most empty Mnl and 
e.g. Mn2,3 atoms (see Table 1) explains the appearance 
of two hole peaks. There is no contribution to the first 
low-energy hole peak from Mnl e g states. They are the 
most empty and contribute only to the second hole peak. 

One would expect that two holes per unit cell move 
apart from each other as far as possible to minimize the 
Coulomb repulsive energy. However the situation is more 
complicated due to the strong p-d hybridization. With 
the help of model hamiltonian describing in the simplest 
way the interaction t between oxygen p level with d levels 
of two neighboring manganese ions 

/ di t\ 
H= d 2 t \ 
\ t t p j 

it is easy to simulate the process of hole distribution 
among the Mn sites. Suppose that one hole is on the 
first Mn site. Then the presence of the hole leads to the 
increasing of the d\ energy due to its smaller occupancy 
and hence to the increasing of the d 2 level energy due to 
hybridization via p oxygen. Then the second hole will oc- 
cupy the g?2 level because it has higher energy. If one hole 
locates on one specific site (let it be Mnl in Fig. || (iii)), 
that leads to decreased electron content on the neighbor- 
ing sites in a (Mn2) and c (two Mn3) directions as it was 
demonstrated for the doped metallic solution. The sec- 
ond hole occupies Mn2 site with higher energy level and 
in its turn makes more empty its nearest neighbors along 
c direction (two Mn4). The hole-rich sites Mnl and Mn2 
are not the same due to essentially different orientation of 
e g orbitals of neighboring Mn atoms along c axis. Such 
system consisting of six hole-doped Mn atoms (Mnl, 
Mn2, two Mn3, and two Mn4) per 2ax26x4c supercell 
is closed and produces small influence on the atoms of 
next a6-slab (Mn5 and Mn6) and on the Mn atoms in 
the back ac-plane. The break of long-range e g orbital 

TABLE I. Contributions of e 9 symmetry wave functions to 
the occupied majority e 3 orbital of different frontplane Mn 



atoms in the 


case of doped 


insulating La7/g 


Sr 1/8 Mn0 3 . 


atom 




<t>x 2 -y 2 


occupancy 


Mnl 


.044 


+.999 


.852 


Mn3 


.257 


-.966 


.941 


Mn5 


.537 


+.844 


1.00 


Mn2 


.117 


-.993 


.941 


Mn4 


.348 


+.938 


.952 


Mn6 


.657 


-.754 


1.00 
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overlapping along c axis leads to the appearance of the 
insulating state. 

The calculated energy gap of 0.03 eV in the doped insu- 
lator is very small and can be smeared out with elevating 
temperature. This temperature may correspond to Too 
or the O* — > O' transition temperature. 

Present orbital order for one undoped plane looks simi- 
lar to that previously proposed in and thus can be used 
for the explanation of their experimental dependence of 
polarization on azimuthal scan. The fact that holes accu- 
mulate on one doped ac plane but there are no holes on 
the other plane is consistent with the observed superlat- 
tice satellite spots in neutron diffraction. The charge 
difference between hole sites and other sites is actually 
small. Since two holes are localized on six Mn atoms, the 
hole charge is 1/3 in average. It was demonstrated earlier 
that hole charge is also spread over oxygen atoms neigh- 
boring the stripe, causing the further reduction of charge 
difference between doped and undoped sites. This may 
be the reason why several experiments could not detect 
the present pattern of charge order. 

It was demonstrated theoretically for a wide range 
of magnetic compounds containing Jahn- Teller ions |l3j | 
that the exchange mechanism (caused by interplay be- 
tween on-site Coulomb interaction U and intersite hop- 
ping t) by itself, without taking Jahn- Teller distortions 
into consideration, is able to give a correct picture of the 
orbital order. This conclusion was supported by band 
structure calculations |l(| of undistorted model PrMn03 
compound. The system Lay/gSriygMnOa is unique real 
compound which has an orbital order without Jahn- 
Teller distortion. And the super-exchange mechanism 
overcomes the double exchange mechanism and provides 
the insulating state with orbital order. 

In conclusions, we have investigated the electronic 
structure of ferromagnetic insulating La7/8Sr 1 / 8 Mn03 
using the LSDA+U formalism. We obtained alternative 
sequence of undoped and doped ac planes with alignment 
of holes in 3-line stripe along a axis. The long period 
orbital order and holes in stripe are consistent with ex- 
perimentally observed properties of this compound. 
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